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ABSTRACT 

The Feinberg-Galanin method for heterogeneous reactors is for- 
mulated by using a two-group model rather than an age kernel. This 
treatment is then extended to take into account secondary effects, such 
as fast fission and thermalization of neutrons inside a rod which may 
contain moderator. The use of a single coefficient in a Feinberg- 
Galanin approach allows one to relate the source and sink strength of 
the fuel element to the thermal flux only. By defining a set of four 
coefficients, it is possible to connect the strengths of thermal- and 
fast-neutron sources and sinks to both thermal and fast fluxes. A 
method is presented for calculation of these four coefficients, a u 
a 2 , and /3 2 . 


I. INTRODUCTION 


The theory of heterogeneous reactors was first devel- 
oped with the homogeneous model by calculating the 
criticality of equivalent homogeneous reactors. The actual 
reactor is divided into cells whose shape is determined 
by the fuel-element distribution. Using the Wigner-Seitz 
unit-cell model, each cell is replaced by an equivalent 
(usually cylindrical) cell whose nuclear parameters can 
be calculated. The heterogeneous cells are then replaced 
in the reactor model by homogeneous regions having 
these calculated nuclear characteristics. 


Subsequently, Feinberg (Ref. 1) and Galanin (Ref. 2) 
developed a heterogeneous method applicable to infinite- 
moderator media. This method was extended to finite 
media of rectangular shape by Meetz (Ref. 3) and to 
those of cylindrical shape by Jonsson (Ref. 4). 

Feinberg and Galanin consider each fuel rod as a 
singularity. These singularities are considered as external 
to the moderator and are treated as localized sources 
and sinks. The properties of the rods are contained in a 
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single constant y, which relates the thermal net current 
to the thermal flux at the surface of the fuel element; this 
allows one to relate the strength of the thermal-neutron 
sink at the rod to the thermal-neutron flux at its surface. 
The singularity is then considered as a sink of thermal 
neutrons and a source of fast neutrons. The strength of 
the source is related to the strength of the sink through 
the coefficient 17, the average number of fission neutrons 
produced per neutron absorbed in the rod. This coeffi- 
cient may be different for each rod. 

For reactors with a small number of solid rods, the 
Feinberg-Galanin method is an improvement over the 
homogeneous method. However, it is not completely 
satisfactory for reactors which may have fuel elements 
containing a significant amount of moderator. In such 
reactors, fast neutrons can be produced in a fuel element 
by fission, but slow neutrons can also be produced by a 
slowing-down inside a fuel element. Neither of the two 
methods noted above treats this kind of reactor properly; 
the homogeneous method does not localize the singu- 
larity, and the Feinberg-Galanin method neglects the 
slowing-down inside the singularity. 

In order to introduce these effects in a heterogeneous- 
reactor calculation, the Feinberg-Galanin method is re- 
constructed in a two-group model (Section II), after 
which the additional source and sink effects in the fast 


and thermal groups are introduced (Section III). Thus, 
each fuel element is represented by: 

( 1 ) a source of fast neutrons ( fission ) 

(2) a source of thermal neutrons ( thermal ization ) 

( 3 ) a sink of fast neutrons ( radiative capture and fast 
fission ) 

( 4 ) a sink of thermal neutrons ( thermal absorption ) 

The source and sink terms are related to each other by 
two constants: ^ and ^ 2 , the average number of neutrons 
produced per fast and thermal absorption, respectively, 
in the fuel element; and by p, the probability that a 
neutron slowing down inside the fuel element reaches 
thermal energy. 

The change introduced here with respect to the 
Feinberg-Galanin method is that the sink terms are re- 
lated to both the thermal and the fast fluxes at the surface 
of the rod. These relationships necessitate the establish- 
ment of four coefficients, which must be determined. 
Section IV outlines a method of obtaining these coeffi- 
cients for a simple fuel element. 

Although the treatment here is for the case of fuel rods 
only, the extension to control rods is straightforward. 
In that case, the coefficients 77 are equal to zero. 


II. FEINBERG-GALANIN METHOD WITH TWO-GROUP MODEL 


Consider an infinite moderating medium containing a 
finite number of fuel elements. These fuel elements are 
assumed to be cylinders of finite length and parallel to 
each other, forming a core embedded in the moderator. 
It is assumed either that the distance between these ele- 
ments is large compared with their transverse dimensions, 
or that the lattice is sufficiently symmetrical that the flux 
near the elements possesses enough symmetry for con- 
sideration of the elements as line sources. In order that 
diffusion theory may be used in the moderator, the 
further assumption is made that the distance between 
two fuel elements is large compared with the diffusion 
length. 


Applying diffusion theory, one can write equations for 
the overall fluxes: 


Tf 

= ?ZS*(p)8(r-r*) (1) 

k = 1 


-D 2 V»* t (p) + 2<‘> * t (p) 

X 

= P2«*,( P ) — ZS*(p)8(r — r*) (2) 
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where <£i(p) and </>o(p) refer to the fast and thermal 
fluxes, respectively, in the moderator; Sr is the removal 
cross section of the moderator; and N is the total num- 
ber of elements. 


On the right-hand sides of Eqs. 1 and 2 are the den- 
sities of the different neutron sources of the system. It 
may be seen that each fuel element is considered as a line 
sink capturing S k thermal neutrons per unit time and 
unit length. The quantity S k depends on the position 
along the z axis of the fuel element, which is the reason 
for introducing S* (p) in Eqs. I and 2. Note that the 
Direc delta functions used in the Equations are two- 
dimensional delta functions, so that the product S k (p) X 
S(r — r k ) is the sink density of the fuel element k at 
point p. 


The geometry of the fuel elements suggests the use of 
cylindrical coordinates; thus, the variables in the differ- 
ential equations are separated. Let 


<if>(p) - 4> (r) <t>(z) 

<f)(z) OC cos B z Z 

R = - - 

2 h 2A z 


(3) 

(4) 

(5) 


where A* is the reflector saving on one side along the z 
axis, along which the reflector need not be infinite. The 
term S k (p), which is directly proportional to the flux, 
can be written as 

S* ( p ) = S* cos B z z 

where S& is now a constant to be determined for each 
rod. Remaining are the two-dimensional equations: 


— V 2 + k\ — -JT X ^( r r *) 

- V 2 ^ 2 (r) + «;*,(r) >(r) 


— LvS*8(r-r*) (7) 

*= i 


where 


1 - Dt 


+ Bl 


K 



(8) 


which is the solution for a superposition of k line sources 
located at different positions r k and satisfies the boundary 
condition for an infinite-moderator medium. 

One can easily find the coefficients A k by making the 
following assumption: The flux at a rod emplacement is 
the sum of a symmetrical rapidly varying function due 
to the rod itself and an unsymmetrical slowly varying 
function due to all the other rods. In computing the de- 
rivative about a rod, the derivative of the slowly varying 
function is neglected when compared with the derivative 
of the other function. The rapidly varying function at the 
rod k is stated as 

<f> k = A k Ko(*i | r - r k |) 

In order to evaluate the coefficients A ky the finite 
radius b of the rod must now be considered. The num- 
ber of fast neutrons leaving the surface of the rod per 
unit length and time is given by 


_ 2 ff foD 1 K,A*K x (<c 1 b) = v S k (10) 

dr 

where S fc is the number of thermal neutrons absorbed 
per unit length and time at the center of rod k; then 


A k 


ySk 

27 rb D i Ki Ki( Kjb) 


( 11 ) 


The thermal constant y, defined in Appendix A, is 
expressed as 

-27rbh(r k ) _ S k 
y fain) <f>*( n) 

where <f> 2 ( r k) and J-tifk) are the thermal flux and net cur- 
rent (taken as positive when directed outward) at the 
surface of fuel element k . The fast flux is thus expressed by 


♦■('> = a,EP,.TK.(«.fc) f> (rt) K " ( '" l ^ r,l) 


(13) 


The thermal equation (Eq. 7) is, then, 

-V s <*, 2 (r) + k\$ 2 (r) = tn£ fo(r k ) Ko(*i |r-r*|) 

fc=l 

8(«‘-r it ) 

k=l 


In the fast-flux equation, the solution sought is of the 
form: where 

f(r)=lA t Ko(« 1 |r-r*|) (9) 


(14) 




rjy_ 


*=i 


m p D 2 2?rf> D l «i Ki(* 1 b) 
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Applying a Fourier transformation to Eq. 14, inverting 
and integrating the result, one finds the following expres- 
sion for the flux: 

fair) = H(|r-r fc |) (15) 

fc— i 


where 

*2 “ *1 

_ (** _ + ~2*Dl) K °( K2 l r_r * : l) ( 16 ) 


III. THE FOUR-COEFFICIENT METHOD 


Considered here is the same assembly as the one 
defined in Section II; now, however, the fuel elements are 
allowed to have non-negligible slo wing-down properties. 
Each fuel element is replaced by: 

( 1 ) a sink of thermal neutrons 

(2) a sink of fast neutrons due to absorption and 
thermalization 

( 3 ) a source of fast neutrons due to fission 

(4) a source of thermal neutrons due to thermalization 
inside the rod 

To write the balance equations, one relates the source 
and sink strengths to the coefficients rj and p. 


4. The Neutron-Balance Equations 

One writes the same equations as those used in Sec- 
tion II; but S^ 1} and S£ 2) , the fast and thermal neutron- 
sink terms, now replace the one thermal-sink term: 

(-D.V* + + 2«>)*(p) =*f)S<»(p)3(r-i*) 

k= l 

+ - !) 2S< 1 >(p) 8(r-r*) (17) 

*= 1 

(-D, v* + 2j 2, )*.(p) = pS»*i(p) 

- £s<’>(p)8(r— r*) +p±S<»(p) 8(r-r») 

(18) 

where 

S£ 1} is the number of neutrons which disappear from 
the fast group per unit length and time at the 
center of rod k, 

S£ 2) is the number of thermal neutrons absorbed per 
unit length and time at the center of rod k . 


rji is the average number of fast neutrons produced 
per fast absorption in rod k ( may depend on k ) . 

t} 2 is the average number of fast neutrons produced 
per thermal absorption in rod k. 

Note the extra terms in the right-hand side of these 
Equations: 

(1) The term (r} 1 — 1) S fc (1) (p) 5 (r— r*) accounts for 
fast absorption and fast fission in rod k . 

(2) The term p S* 1} (p) $(r— r k ) accounts for ther- 
malization in rod k; the overall resonance-escape 
probability p is used as an approximation. 

In these Equations, the z dependency can be removed by 
letting 

*(P) = 4>( r ) cos B z z 

Then, 

(-V* + (r) =-2l£s<»8( r-r*) 

jfc=l 

+ 5l D^E S r , 8(r-r i ) (19) 


(~V« + «■)*(,) = -£-Z* *(r) 

- -T t S l 2> '(r-r>) +£± S<» 8(r-r») 

U 2 k= 1 U * fci 


where 


( 20 ) 


VfXi 

Dx 

2j 2 ’ 


a) 


+B* 


D 2 


+B 2 


B: = 


2 h + 2a, 


( 21 ) 


It is now necessary to express S' 1 ' and S ( h 2) in terms 
of the fluxes at the surface of each rod <f> i ( r* ) and <f> t (r k ). 
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Therefore, the coefficients are defined: 

«! is the probability that a fast neutron entering the 
fuel element escapes from the fuel element as a 
fast neutron. 

j$ x is the probability that a fast neutron born from 
fission inside the fuel element escapes from the 
fuel element as a fast neutron. 

a 2 is the probability that a thermal neutron entering 
the fuel element escapes from the fuel element. 

fi 2 is the probability that a neutron thermalized in- 
side a fuel element escapes from the fuel element. 

Using these coefficients, one can write: 

S<» = 2mbn (1 - “0 + ^ S I fc 2> _ ft) 

+ , 1 S< 1 > (1 -/?,) (22) 

S‘» +ps<«(l -fit) (23) 

where f is the partial current going inward at the rod 
surface, and b is the radius of the rod. 

The fast and thermal net currents can be expressed by 

2ir bj, = 2nrb (/+ - /") = V2 S k 2) ft + S k 1 ’ (^ift — 1 ) 

(24) 

2nb J 2 = 2jrfc ( /' 2 -jl) = P Sl l} ~ Sr (25) 


where 


— , (1 — «i)(l + «z) 


fli — irb 


M 


(l - a 2 )(30! - 2 - cc^O 

b, — — 7rb 7/2 ^ 

i (1 a 2 ) [1 + iyi(3)3i — 2) + ai(l ~ qifti )] 

a 2 — irb - jyf 

„ (1 - 00(3/3, -2- a*/9 a ) 

bz — — irb p ± jj 

M = [1 + 171(3^1 — 2) + a,(l — 171JS1)] (1 + a 2 ) 

— j>» 7 2 (3/?2 2 a 2 /? 2 ) (3/3 i 2 a i^i) 

(30) 

Thus, the differential equations ( Eqs. 19 and 20 ) become 
( - V 2 + ki) ^(r) = X) (<£i*) 8(r - n) (31) 

(- V 2 + ) *,( r) = -g-S E 

+ £(**) S(r - r*) (32) 

where 

(<#>ifc) = -^-{^>i( r *) t^ft + (>?» - 1 ) a >) 


+ U 20-2 + (yi — 1) £*i] ) 


(33) 


(<j>2k) — T^“{$i( r *) [p° 1 ft) <t>2( r k) [pbi a a]} 

(34) 


and, from diffusion theory, 

n + U = ( 26 ) 

(27) 

where the ;’s and </>’s denote the partial currents and fluxes 
at the outside surface of a rod. 

Note that, because of the assumption of symmetry 
made at the beginning of Section II, the flux at the 
surface of a rod is independent of the azimuthal angle. 

With the six linear relationships stated in Eqs. 22-27, 
one can find the S k ’s in terms of the </>’ s: 

Sfc 1 * = fli <£i(r*) + ft <f>i(r k ) (28) 

Sl 2> = a 2 </> 2 (r k ) + b 2 <f>i(r k ) (29) 


B. Solution of the Balance Equations; 
Criticality Condition 

For an infinite moderator, a solution of Eq. 31 is 


*i(r) = E A,K 0 (K 1 |r-r k |) (35) 

k= 1 

To compute A k , one again neglects, in the derivative of 
the flux at a rod k, the derivative of the flux due to the 
other rods when compared with the derivative of the 
flux due to the rod k. The current at the boundary of the 
k element is computed in the same manner as that used 
in Section II, with the result: 


4k = 


where 


(<ftifc) 

2>rfe K 1 K l (K 1 b) 

1 


(</>ifc) F 


F = 


2 irb k 1 Ki(kJ}) 


(36) 

(37) 


5 



JPL TECHNICAL REPORT NO. 32-333. 


Then, 

fc(r) =Fi:(* 1 *)Xo(« 1 |r-i*|) (38) 


By substitution into Eq. 32, 


(- V* + K 2 2 )<^ 2 (r) k-r|)+ £ (**) 8(r - r*) (39) 

As in Section II, a Fourier-transform procedure is used to solve this Equation, resulting in 

<h( r ) = ~D~—; i" fiF 2 E (^) [Ko(/fi |r - r*|) - K„ (« 2 |r - r*|)] 
u 2 k 2 — Kl 

+ "2^ E(4>2*) K 0 (k 2 |r — r*|) (40) 

k - 1 

To obtain the criticality condition, one computes from Eqs. 38 and 39 the thermal and fast fluxes 
at the surface of each rod; 

<h(r m ) — -jL T [<f> i(r*) [T] 2 b 2 + ( Vl — l)aj + <£ 2 (r t ) — l)b,]} K 0 (*i [r M — r fc |) 

( 41 ) 

<^ 2 ( {£»(**) U2&2 + (>;i — l)aj 

+ Qzifk) [v&i + (iji — l)hj]} (Ko (*i|r m - r*|) - K 0 (k 2 |r m - r*J)} 

+ "^T {</*i( r *) Cpo i - W + *,(ifc) [phi - a 2 ] } K„(k 2 |r m - r*|) (42) 

The result is a system of 2N linear homogeneous equations whose 2N unknowns are the <j>i(r k ) 
and <f> 2 (r k ). Non-trivial solutions exist only if A = 0 (or, the determinant of the system is zero), The 
criticality condition is that the 2N-order determinant A must vanish. 

In practice, use can be made of symmetrically situated fuel elements having the same surface 
flux <t>(r k ), and the number of unknowns can thus be considerably reduced. Most often, the order 
of the determinant is smaller than 2 N, 
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IV. CALCULATION OF THE FOUR COEFFICIENTS 


The four coefficients a^ou/^, introduced in Section 
III, depend on the nuclear properties and geometrical 
configuration of the fuel element. Each fuel element, if 
different from the others, can have a different set of these 
coefficients. Considered here is a simple model of a fuel 
element, accompanied by a method of obtaining these 
coefficients. 

The element is composed of two concentric cylindrical 
regions. The inner region contains moderating material 
(2 a < 2,), and the outer region contains fuel (2 a ~ 2,). 
The cross section of the fuel element is shown in Fig. L 



Fig. 1. Cross section of fuel rod containing moderator 


First, a few probabilities are defined which will be 
useful in obtaining these coefficients: 

?! is the probability that a neutron coming from 
the inner moderator goes through the fuel shell 
without making any collision. 

P 2 is the probability that a neutron escapes from 
the fuel element after a scattering collision in 
the fuel. 

P 3 is the probability that a neutron enters the inner 
moderator after a scattering collision in the fuel. 


7 ri is the probability that a neutron coming from 
outside goes through the fuel shell only and 
escapes from the fuel element without a collision. 

tt 2 is the probability that a neutron coming from 
outside reaches the inner moderator without a 
collision in the fuel. 

t r 3 is the probability that a fast neutron (born in 
the fuel annulus) entering the inner moderator 
does not thermalize inside the moderator. 

tt' is the probability that a fast neutron (born out- 
side the fuel annulus) entering the moderator 
does not thermalize inside the moderator. 

t r 4 is the probability that a thermal neutron escapes 
the inner moderator. 

These probabilities (except 7r 3 , and 7 r 4 ) must be de- 

fined for both fast and thermal neutrons. In general, 
accurate calculation of these transmission probabilities 
is quite difficult and could be the subject of a separate 
study. For the present purposes, it was considered suf- 
ficient to make rough calculations utilizing rather gross 
simplifying assumptions. These assumptions and the cal- 
culations are outlined in Appendix B. The probabilities 
are defined graphically in Fig. 2. 

The following analysis shows how the detailed trans- 
mission probabilities are combined to give the four coeffi- 
cients desired. 


A. The Coefficient a,: Probability That a Fast 
Neutron Entering the Fuel Element 
Escapes as a Fast Neutron 

Referring to Fig. 2, assume that S fast neutrons enter 
the fuel element per unit time and per unit length, and 
that x of these neutrons, after entering the inner modera- 
tor, leave it again as fast neutrons. Hence, 

[S ( 1 ir i -*>) + X(1 - P i)l « 

neutrons make at least one scattering collision in the fuel, 
whereas 

[S(l 7T 2 ) + X(1 - Pi)] « (1 - P* - P 3) 
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Fig. 2. Escape and transmission probabilities 
in fuel model 


make two scattering collisions in the fuel, and so on. 
Here, w — The total number of scattering colli- 

sions these neutrons make in the fuel is (by summation) 


[S(l — 7Ti — 7 r 2 ) + x(l — Pi)] (X) 

1- (l-p,-p.)« 


Thus, the number of escaping neutrons is 


Sax = Sir, + £Hl^-=_Z-4 + ^ 1 ~ P ^ ]m P, + xFl 

L ~ (a — r 2 — j r 3 j w 

(43) 


and the number of neutrons entering the inner moderator 
is 


where 




7T 3 


7T2 (1 _ ?2 ~ P. 3)^2 + ( 1 7T 1 7r 2 )wF 3 

1 - (1 ~ P 2 - P 3 )<o - (1 - PJ^Pair; 


(46) 


8. 77ie Coefficient p t : Probability That a Fast 
Neutron Born From Fission Inside the Fuel 
Element Escapes From the Fuel Element 
as a Fast Neutron 


Assume that S neutrons are born from fission in the 
fuel element per unit length and time, and that x of 
these, after entering the inner moderator, leave it again 
as fast neutrons. Then, 

[S(l -P 2 ~ P 3 ) Hb x(l - P x ) ] w 

make at least one scattering collision in the fuel. The 
total number of scattering collisions these neutrons make 
in the fuel is 


[S(l - P 2 - P 3 ) + x(l -P,)]« 

1 - (1-P a -P B )« 

Thus, 

Sfr = SP 2 + [ - (1 7 + **>. 

(47) 

of these neutrons escape, whereas 

x „„ . [S(l - P, - P 3 ) + x(l - P,)]< 

1 - (1 - P 2 - P 3 )<* 


— = SP 3 + 

^3 


-P» 

(48) 


neutrons enter the inner moderator. Hence, 
(1 P2 P 3 )<oP 2 


181 = P 2 + 

+ 


1- (1-P,-P,)« 

P 3^3 


X 


1 (1 P-2 P 3 )«) (1 Pi)(l>P 3 W 3 

P, +<o(P 2 + P,P 3 - Pi) 

1 - (1 - p 2 - p,)» 


(49) 


-4 = S,r 2 + 


[3(1 ir, 7r 2 ) + x(l — P, ) ] 1 
1 -(1 -P 2 -P 3 )«, 


(44) 


Solving the two previous Equations for a, gives 

_ (1~ IT, — 7T 2 )<»P 2 

1 ^ 1 - (1 - P 2 -P,)» 

(P 2 + P,P 3 ~ Pl)(» + P, 

1 1 -(!-?,-?,)• 


(45) 


C. The Coefficient a,; Probability That a Thermal 
Neutron Entering the Fuel Element 
Escapes From the Fuel Element 


One can use here the same procedure as for a,, where 
vt replaces ir' . This yields 


CU — 7Ti + 


+ *2 


(1 7Ti 7T 2 )<dP 2 

1 - (1 ~ P 2 - P 3 ) 01 

(P 2 + P 1 P 3 -P 1 ) tt , + P, 
1 - (1 -p,- P 3 )cu 


(50) 
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where 


e 2 — 7 r 4 


7T 2 (1 — P 2 _ P 3 )<^2 + (1 

1 - (1 “ P. " Pa)® ~ (1 


7Ti 7T 2 )t0p3 
P 1 )<t>P 3 7T 4 


( 51 ) 


D. The Coefficient p 2 : Probability That a Thermal 
Neutron Thermalized in the Fuel Element 
Escapes From the Fuel Element 

Assume that S neutrons are thermalized in the element 
per unit length and time, and that x of these, after enter- 
ing the fuel, return into the inner moderator and sub- 
sequently escape again. Then, (S + x ) (1 — neu- 
trons make at least one scattering collision in the fuel. 
The total number of scattering collisions these neutrons 


make in the fuel is 


(S + x) (l-Pi)o> 

1 - (1 - P 2 - P 3 )c 


Thus, 


_ /c _i - r \ p -u (S + x) (1 Pi)o> p 
S /?2 - (S + x) Pi + i _ (1 _ p 2 — P 3 )a/ 2 


( 52 ) 


neutrons escape from the rod, and 

x (S + xHl-PQo) p 
tt 4 1 - (1 - p 2 - p*h 3 

neutrons re-enter the inner moderator. The result is 

Pi + 0>(P 2 + P1P3 - Pi) 

^ 1 - (1 - P 2 - P 8 )« - (1 - Pl)wP 3 7T 4 


( 53 ) 


( 54 ) 
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V. CONCLUSIONS 


This development of the heterogeneous method, using 
the four-coefficient and two-group approach, extends the 
area of application of the method to include a greater 
variety of elements, with no expected decrease in the 
order of accuracy. The same characteristics are present 
as in the single-coefficient technique; the increased com- 
plexity is achieved at the cost of increased computation 
time. 

The size of the criticality determinant is twice that of 
the same determinant in the single-coefficient case. It is 
therefore advisable, before using this technique on a 
complex problem, to examine the array for all possible 
symmetries. In a regular hexagonal lattice, for example, 
reduction of the magnitude of the determinant by a fac- 
tor of 12 is possible. 

As presented, the calculation of the individual trans- 
mission probabilities is a separate problem and can be 
performed to any order of precision and complexity re- 
quired. Multigroup-cell theory or Monte Carlo methods 
could be utilized if necessary, rather than the simple 
model presented here. 

The flux distribution in the moderator can be obtained 
from Eqs. 38 and 40. These Equations include the terms 
and ( (j> 2 k)y defined in Eqs, 33 and 34. They depend 
on <£i(r fc ) and <j> 2 (r k ) } the fluxes at the surface of each 
fuel element. The latter are solutions of the homogeneous 
linear equations (Eqs. 41 and 42), whose determinant has 
been made zero. Hence, the solution is not unique. A 
value for one <£(r*) must be arbitrarily chosen. By sub- 
stituting this value in Eqs. 41 and 42, one gets a system 
of linear non-homogeneous equations. One of these equa- 
tions is redundant, because one more equation than un- 
knowns remains. This non-homogeneous system can be 
solved by calculating determinants of order 2 N — I, at 
most. 


A program has been drawn up for the IBM 7090 com- 
puter which will accomplish both calculations of the 
coefficients and solution of the determinants. 

Further improvements can be made in the four-coeffi- 
cient method. As was pointed out above, the resonances 
in the fuel are taken into account by an overall resonance- 
escape probability. Actually, since the resonance absorp- 
tions take place mainly in the fuel, it might be worthwhile 
to take into account the localization of this phenomenon. 
This might be done, for instance, by introducing a third 
group. One observes that the resonance absorptions occur 
primarily in a relatively narrow range of energy for 
uranium (between about 2 and 200 ev). Using the flux 
distribution of this intermediate group, it is possible to 
determine at what rate these neutrons flow into each 
fuel element; also, it is possible to evolve a set of 
cross sections for this group which includes resonance 
properties. 

The introduction of three groups requires the use of 
nine coefficients, rather than the four employed in the 
method outlined in Section III. These nine coefficients 
relate the three sink terms to the three fluxes at the 
surface of each fuel element and are equivalent to the 
four coefficients defined in Eqs, 28 and 29. 

The next improvement would be to consider a finite 
medium, which is necessary for the case of thin reflectors. 
For example, one might apply the refinements developed 
in the present investigation to the Jonsson theory ( Ref. 4 ) . 

It should not be forgotten that all these extra effects, 
if taken into account, give rise to mathematical compli- 
cations. For instance, the size of the determinant of the 
criticality condition is tripled if three groups are used. 
Hence, depending on the number of fuel elements to be 
used, a compromise must be found between the degree 
of complexity one can afford and the accuracy one 
desires. 
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APPENDIX A 


The Thermal Constant in Cylindrical Geometry 


The thermal constant 7 of a slug in a diffusing medium 
is defined as the ratio of the net flow per unit time of 
thermal neutrons into the slug to the value of the thermal 
flux at the surface of the slug. 


One way of obtaining an approximate value of this 
constant is presented here, based on the following 
assumptions: 

( 1 ) The fuel rod consists of only one kind of material. 

(2) The angular distribution of neutrons entering the 
slug is isotropic (a first-order correction to this 
approximation is performed here). 

(3) The collision density is constant inside the fuel 


The thermal coefficient is related to the transmission 
coefficient r, or the fraction of all neutrons incident upon 
the surface of a lump which pass through the lump with- 
out being absorbed: 


T 


y — 


L 

f 

- 2wh] 

4> 


= *b ] - 


r 


— 7T b 


r + r 

1 T 

1 H~ T 


( A-l) 


The terms J and <f> are the net current and the flux, 
respectively, at the surface of the rod; j + and j~ are the 
partial currents at the same surface ( the f current going 
outward ); b is the radius of the slug, which is assumed 
to be circular cylindrical. From this relationship, t is 
computed, and A is deduced. 


1. The Transmission Coefficient r 

Among S neutrons entering the slug, some scatter, 
some are absorbed, and some escape without undergoing 
collisions. In paragraph 2 of this Appendix, a computa- 
tion is made to determine AS, the number of neutrons 
which escape without undergoing any collision. The 
number of neutrons (1 - A) S, therefore, make at least 
one collision each. Let 

2 t = the total macroscopic cross section of the material 
inside the lump 


2, — the macroscopic scattering cross section of the 
same material 

Thus, among the (1 — A) S neutrons which make 
collisions, (1 - A) (2*/2 r )S make scattering collisions. 
In paragraph 3 of this Appendix, a computation method is 
given for £, the probability that those neutrons which 
have made scattering collisions then escape from the slug. 
For that calculation, the two following assumptions are 
made: 


( 1 ) In the laboratory system, the scattering is isotropic 
in the slug material. This assumption is acceptable 
for heavy nuclei where JL 0 = 2/3A. The average 
cosine of the scattering angle becomes small when 
compared to unity. 

(2) The collision density in the slug is independent 
of position. 

Once £ is found, it is seen that 

(l-x)|s-(l-C)S 

of the incoming neutrons make at least two collisions 
each, and that 

(1-»)£<1-C>%s 

neutrons make at least two scattering collisions each. 
By the same procedure, it is evident that 

neutrons make at least ( n + 1 ) scattering collisions each, 
and that, at each scattering generation, a fraction £ of the 
scattering neutrons escapes from the slug. 


Thus, the total number of neutrons which escape from 
the slug after any number of collisions is expressed by 

St = L+ (1-A)|jc+ (1 -A)^(l 


(A-2) 
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Hence, from Eq. A-l, the thermal constant is 


y 


7 (3- ~ a) (S?. — S«) 

(1 + A) (2 r -2.) +2& s 


( A-3) 


2. The Probability X 

The term A is the probability that a neutron entering 
the slug will escape without making any collision. This 
probability is equivalent to the transmission coefficient r 
of a lump, computed in a first-flight approximation ( Ref. 5, 
p. 247); however, the total cross section is used rather 
than the absorption cross section. 


a. Expression for the probability X. Needed here is the 
angular distribution of the velocities of the thermal neu- 
trons impinging upon the surface of the slug (Fig. A-l). 
A first-order correction to the isotropic distribution is 
given by the diffusion theory (Ref. 5, p. 171): 

/-(/*,*) dftdt =£^[*(0) + (A-4) 

where 

fx is the cosine of the angle which the direction of 
motion of the incident neutron makes with the 
normal to the surface of the lump. 

<\> is the azimuth of the neutron direction about the 
normal. 


2* w) is the scattering cross section of the diffusion 
medium outside the lump. 

p is the angle between the direction of motion of 
the incident neutron and the flux gradient. 

In the present case, the flux gradient can be considered 
as always oriented perpendicular to the surface of the 
slug, neglecting the gross variation of the flux in the 
reactor when compared with the local variation at the 
slug boundary. 

Hence, cos P — fi, and the neutron angular distribution is 

fW) d/id<l> = -[+(0) + (A-5) 

where V^(0) can be approximated from a diffusion- 
theory calculation. 


For cylindrical geometry, 

V4>(0) _*'(<>)_. *h(*b) 

*<0) *(0)- I 0 (*&) 


(A-6) 


where the slug-material constants are expressed by 


a 

D iF) 


(A-7) 



Fig. A-l. Projected path of neutron impinging on 
surface of fuel slug 


Hence, 

fit 1 * 4>) d(xd<t> — <£ 0 (1 -|- g^) (A-8) 


where 


K Ii(xb) 

g ~Y^ljTb) 


(A-9) 


The integral, over all directions into the lump, of the 
velocity distribution is 


S'/. 1 ^ + g ^ dii /" ^ = ^r( x + ¥) (A ' 10) 

The probability that a neutron will pass through a distance 
of material without making any collision is e~^s y where 
is the total cross section of the material, and s is the 
path length shown in Fig. A-l. 

The probability A can now be expressed as 


A = 


f f 

J J /(p» e dfid4 


j f J f i(p> 


(A-ii) 


</>) djxd<f) 
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A = 


I 


( 


-f) 


f ^(1 + gp) dp I # 


2x 

e d<f> 


(A-12) 


b. Calculation of the probability X in cylindrical geom- 
etry from Eq, A-12 . Referring to Fig. A-l, one makes 
the change of variables defined by 


= cos ft cos y 


tan = 


tan y 
sin ft 


s = 


2b 


cos ft 

COS y 


The Jacobian of the transformation is 


J = cosy 


(A-13) 


Hence, 


A = 


( i+ f g ) 


X 


( f r */2 „ , "2r2fc( cos /S/cosy) 

<JJ cos y3 cos 2 ye a fid, y 


+ g 


/r 


cos 2 cos 3 y e 


-2 t2£>(cos £/cos 7) 


} 


dftd y 

(A-14) 

These integrals were evaluated numerically. Writing, 


2 fb 



Fig. A-2. Numerical values of f s and l 2 vs fuel-rod 
radius measured in mean free paths 


A = 



(i i + gh) 


the values of I x and I 2 are plotted in Fig. A-2 as functions 
of H T b. 


der, and £ is the probability that a neutron bom in this 
source will escape. 

This problem has been treated by Cohen and Estabrook 
(Ref. 6), with the following result: 


3. The Probability £ 

The term £ is the probability that a neutron will escape 
from the cylindrical slug after a scattering collision. It is 
assumed that the scattering is homogeneous and isotropic 
inside the slug. Hence, the system is equivalent to a 
uniform-source material whose shape is a circular cylin- 


f = ^ j- 2 + ( 2 + jL) Jx(2^) K.i'Zrb) 

+ I 0 (2 r fc) Kx(2 r fe) -hCZrb) Ko(2 r &) 
+ 22 T fc/ 0 (2 r fc)Ko(2 r i)( (A-15) 
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APPENDIX B 

Approximate Method for Evaluating Neutron-Transmission Probabilities 


In computing the transmission probabilities, the follow- 
ing approximations are made: 

( 1 ) The scattering is isotropic in the laboratory system. 

(2) The angular distribution of neutrons impinging on 
the outside surface of the fuel element is isotropic. 

( 3 ) The collision density 2 r </> in each region of the fuel 
element is space-independent for both thermal and 
fast neutrons. 

1. The Probability P x 

The term P x designates the probability that a neutron 
coming from the inner moderator goes through the fuel 
shell without making any collision. It is assumed, as a first 
approximation, that the angular distribution of the neu- 
trons entering the moderator is uniform. The probability 
that a neutron coming out of the element of area dA 
(Fig. B-l) through a solid angle dil about 0 and <f> goes 
through the fuel shell without colliding is e'ZrSo where 2 r 
is the total cross section of the fuel. With use of the 
approximation above, it follows that the expression 

dn = —cos 2 0 cos cj> d<f> d$ 

IT 

is the probability that a neutron going through dA is in 
a solid angle (kl about ( 9 and <f> (1/tt is a normalizing 
factor). The product e^rSo dQ is the probability that 
a neutron going through dA is in a solid angle dn and 
penetrates the fuel shell without a collision. The proba- 
bility P l? the sum of these elementary probabilities over 
the solid angle, is given by 

Pi = — n /2 cos 2 0 cos <f> e~^ TSo cLpdO (B-l) 



Fig. B-l. Three-dimensional path of neutron from inner 
moderator passing through fuel shell 
without collision 


where, according to Fig. B-l, 

R 

cos 6 


So 


One may then write 

P, 


('V l “F sinI *~I cos *) 


2. The Probability P 2 

The probability that a neutron escapes from the fuel 
element after a scattering collision in the fuel is denoted 
by the term P 2 . On the assumption that scattering is 
isotropic in the fuel, and that the flux can be considered 
(B-2) as a constant inside the fuel, each volume element dv is 


where 



f v/2 r-tr /2 

Jo Jo 


cos 2 6 cos <£ e 


-2 t (R/ cos 9 ) ^ 1 - Bin 2 £--~coa 


d<j> d$ (B-3) 
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considered as a unit source. The probability that a neu- 
tron coming from dv reaches an area dA on the outside 
surface of the fuel element is the product of e^rs and 
the fraction of solid angle through which dA is seen from 
dv, or 

dA cos 0 cos </> 

4 ttS 2 

where s is the distance between dA and dv. 

To obtain P 2 , this product is integrated over the volume 
and divided by the source strength, which is v: 

P 2 = — f f f f — C ° 3 * S ^ COS — e s 2 ds cos 0d4>dO 
VJa j* j e ArrS 2 

(B-4) 


Note that there are two kinds of limits of integration, 
depending on whether <j> is larger or smaller than a (Fig. 
B-2) where a — sin -1 (r/R). These two limits are 




Therefore, the resulting expression is 


Fig. B-2. Three-dimensional paths of neutrons escaping 
from fuel through outer wall 



where 


3. The Probability P 3 

The term P 3 represents the probability that a neutron 
enters the inner moderator after a scattering collision in 

the fuel. Considered now is a small volume element dv 
in the fuel and a small area dA on the surface of the 
inner moderator. Applying the same reasoning as for P 2 , 
the limit of integration ( Fig. B-3 ) is now 



w (4 . R , = r /2 I ° COS 2 e cos 4, e ' 2r<B/c08 ” [ C08 * " V(^) “ - 8in ’ 84 dfi 
\R / 0 *^=0 


15 




JPL TECHNICAL REPORT NO. 32-333. 



Fig. B-3. Three-dimensional path of neutron escaping 
from fuel to inner moderator 

and the result of the integration is 



where z ( r/R, R, 2 T ) is defined by Eq. B-3. 

4. The Probability i T t 

The probability that a neutron coming from outside 
goes through the fuel shell only and escapes from the fuel 
element without a collision is represented by 7r t . Assume 
an isotropic angular distribution for the velocities of the 
neutrons coming into the fuel element, 

/ '( 0, <£ ) d<f> dO — cos 2 0 cos <p d<f> do ( B-9) 


The probability that a neutron traveling in a solid angle 
dQ about 0 and <f> escapes is provided that <f> > a ( see 
Fig. B-2). Hence, 

_ // </>) e^^d4_de 
j “( 0 , <p ) d<j> dd 

Since 

s 2 — 2R 

COS 0 

then, 

Vi = ~ I (i ,R,2r ) (B ' 10) 

where I (r/R, R, 2 r ) is defined by Eq. B- 6 . 

5. The Probability 7r 3 

The term 7 r 2 denotes the probability that a neutron 
coming from outside reaches the inner moderator without 
a collision in the fuel. Using the same procedure as for 
ttj, one obtains 

= e ~^_d4_de 

f / j~(6, 4>) d<f> d8 

where, according to Fig. B-2, 

Therefore, 

T ^T W (l’ R ’ Sr ) (B ' n) 

where W (r/R, R, 2 T ) is defined by Eq. B-7. 

6. The Probability 7 r 3 

The term tt 3 expresses the probability that a fast neu- 
tron entering the inner moderator does not thermalize 
inside the moderator. Two cases are distinguished: 

(1) The neutrons born inside the slug from fission. 
These neutrons have a well known average leth- 
argy, and one can have an idea of the average 
number of scattering collisions which will make 
these neutrons thermal. 

(2) The fast neutrons which enter the slug from out- 
side. These neutrons belong to the fast group, but 
actually their lethargy is not well defined and is 
spread between thermal and minimum lethargies. 
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In case (1), knowing the probability that a neutron 
escapes at each scattering collision, one can determine 
the probability that a fast neutron escapes without be- 
coming thermal. From the notations and results of Ap- 
pendix A, 

In these expressions, 


A + (1 - A) (l 



t 


(B-13) 


A = the probability that a neutron going through the 
moderator does not collide. 


A = -/(0.f,S r ) 

7 r 


£ = the probability that a neutron, after a scattering 
collision in the moderator, escapes from the mod- 
erator. 


£ = - 2 + + t) (*> (*) + *0 (*) Ki (*) 


In this case, it can be assumed that 2 a /2 r — 1, because 
only fast neutrons are considered. 


(I) Ko (I) 4- 22I 0 (I) Ko (I) 


Calling N aV g the average number of collisions which 
make a fission neutron thermal, one obtains 

irs = A+ (1- A)C+ (1-A) (1-0C + 

"’ + (1“ A)(l” £) y ar? 

_ A i /i — y \ ~ H ~ £)* gr g 
71,3 A +■ ( 1 A j £ 1— (1 £) 

= A+(1-A) P-U 

(B-12) 

Turning to case (2), in the context of the two-group 
model, one can say that at each scattering collision the 
average probability that a fast neutron becomes thermal 
is 2^/2*, neglecting fast absorptions in the moderator. 
Therefore, using the results of Appendix A, 


where 

l~ r 

7. The Probability tt 4 

The probability that a thermal neutron escapes the 
inner moderator is represented by 7r 4 . This probability 
accounts for the thermal absorption which takes place in 
the moderator. Most of the time, it is very close to 1. The 
term 7r 4 is equivalent to the transmission coefficient of 
the moderator slug. Hence, according to the results of 
Appendix A, 

W * = X + (1-a)^-— — ( B _ 14 ) 
where 2 a /2 r of the moderator is used. 


NOMENCLATURE 


b radius of a fuel element 
B z longitudinal buckling 
D u D 2 fast and thermal diffusion constants 
h half-length of a fuel element 

K n (x) modified Bessel function of the second 
kind and of the nth order 

N total number of fuel elements 

Ni average number of collisions to make 
a fission neutron thermal 


p resonance-escape probability 

r radius of the inner moderator cylinder 
of a fission-electric cell element 

r two-dimensional vector: space variable 
in transverse flux equations after 
removal of the z dependency 

R outside radius of the fuel layer of a 
fission-electric cell element 

S< 1} , S£« fast and thermal sink strengths of the 
singularity 
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NOMENCLATURE (Cont’d) 


a sin -1 ( r/R ) 

a,, ft 2 four coefficients connecting sink 

strengths to fluxes 

y thermal constant 

reflector savings on one end along the 
z axis 

€ scale factor 


> ^a 2) f as * an< 3 thermal macroscopic 
absorption cross sections of the 
pertinent material 

2^ macroscopic removal cross section 
of the moderator 

2* macroscopic scattering cross section of 
the pertinent material 


£ probability that a neutron, after 
colliding in the fuel shell, escapes 
without making an extra collision 

■q x average number of neutrons produced 
per fast absorption in the fuel element 

average number of neutrons produced 
per thermal absorption in the fuel 
element 

#c transverse component of the inverse of 
diffusion length 

A probability that a neutron going through 
the fuel shell does not collide 


2 r 

<M r ) 

Mr) 

<t>i( r k),M r *) 


macroscopic total cross section of the 
pertinent material 

fast flux at point r 

thermal flux at point r 

flux at the surface of the fuel element 


(p) three-dimensional vector: space 
variable in overall flux equations 

t transmission coefficient of a slug 

o> ratio of the scattering and total cross 
sections in the pertinent material 


REFERENCES 


1. Feinberg, S. M., “Heterogeneous Methods for Calculating Reactors/* Proceedings 
of the International Conference on the Peaceful Uses of Atomic Energy, Geneva, 
Switzerland, Vol. 5, Paper No. 669, p, 484, 1955. 

2. Galanin, A. D., “Thermal Coefficients in Heterogeneous Reactors,” Proceedings of 
the International Conference on the Peaceful Uses of Atomic Energy, Geneva, 
Switzerland, Vol. 5, Paper No. 666, p. 477, 1955. 

3. Meetz, K., “Exact Treatment of Heterogeneous Reactors,” Proceedings of the Sec- 
ond United Nations International Conference on the Peaceful Uses of Atomic Energy , 
Geneva, Switzerland, Vol. 16, Paper No. 968, p. 61 1, 1958. 

4. Jonsson, A., and Naslund, G., “Heterogeneous Two Group Diffusion Theory for a 
Finite Cylindrical Reactor,” Aktiebo/ogef Atomenergi AE-57, Stockholm, 1961. 

5. Meghreblian, R. V., and Holmes, D. K., Reactor Analysis , McGraw Hill Book Com- 
pany, Inc., New York, 1960. 

6. Cohen, E. R., and Estabrook, F. B., Surface Fluxes and Currents for Various Shielded 
Radiation Sources, AECU-3595, Atomics International Div., North American Avia- 
tion, Inc., Canoga Park, Calif., 1957. 


18 



